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• Premise of the study: Microsatellite markers were developed and cross-species transferability assessed for two Persoonia spe- 
cies to evaluate genetic diversity and population genetic structure of these broadly distributed southwest Australian tree 
species. 

• Methods and ResuUs: Microsatellite-enriched libraries and 454 GS-FLX shotgun sequencing were used to identity nine micro- 
satelhte loci for P. elliptica (one 454; eight cloning) and six for P. longifolia (three 454; three cloning). These loci were 
screened for variation in individuals from populations in southwestern Australia. In P. elliptica, observed and expected 
heterozygosities ranged from 0.46 to 0.93 and 0.42 to 0.88, respectively. For P. longifolia, observed and expected heterozygosi- 
ties ranged from 0.04 to 0.88 and 0.04 to 0.84, respectively. 

• Conclusions: The microsatellites identified in this study will enable the examination of population and spatial structuring of 
genetic diversity in P. elliptica and P. longifolia, two priority species for mine site restoration in southwestern Australia. 
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Persoonia Sm. (Proteaceae) is an endemic Australian genus 
of woody perennial plants containing 100 species, 43 of which 
are found in Western Australia (Weston, 1995). Persoonia el- 
liptica R. Br. and P. longifolia R. Br. are widespread, small trees 
within the southwestern Australian jarrah forest and are priority 
species for mine site restoration practitioners. These two sym- 
patric congeners are key components of the jarrah forest and 
share pollinators and seed dispersers, but contrast markedly in 
their population densities. There are no species- specific molec- 
ular markers available for these two diploid species. Here we 
report the isolation and characterization of microsatellite mark- 
ers for both species using (1) cloning of microsatellite-enriched 
libraries {P. elliptica) and (2) 454 GS-FLX shotgun sequencing 
{P. longifolia). This will enable the examination of genetic di- 
versity, range-wide genetic differentiation, and mating system 
parameters in the two species. 

METHODS AND RESULTS 

Fresh leaf material was collected from P. elliptica plants in the southwest 
of Western Australia (population code: Pe-AE, 31.572°S 116.181°E, collec- 
tion number: JS102; population code: Pe-AW, 31.613°S 116.151°E, collection 
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number: JS102). Herbarium material is deposited at the Kings Park and Botanic 
Garden Herbarium (KPBG), Perth, Western Australia, Australia. Fresh leaf ma- 
terial was also collected from P. longifolia plants in southwestern Western Aus- 
tralia (population code: Pl-MD, 32.176°S 116.245°E, voucher: R. Davis 3717 
[PERTH 0481230]; and population code: Pl-CE, 33.356°S 1 16.345°E, voucher: 
A. Gundry 486b [PERTH 04801881]). The fresh leaf material was stored at -80°C 
and DNA extracted from the frozen material using the procedure described 
in Jobes et al. (1995) with the following modifications: 4 |aL of proteinase-K 
was added prior to placing samples in the water bath and DNA was precipitated 
in the final stages by added 95% ethanol rather than isopropanol. 

Genetic Identification Services (La Canada Flintridge, California, USA; 
http://www.genetic-id-services.com/) was employed to develop microsatellite- 
enriched libraries for P. elliptica for four different repeat motifs (CA^, GA^, 
ACCn, and ATGJ. Briefly, genomic DNA was restricted with seven blunt-end 
cutting enzymes {Rsal,HaeUi, BsrQl, PvwII, Stul, Seal, £^coRV). Fragments in 
the size range of 300-750 bp were linker adapted with oligonucleotides that 
contained a Hindlll site and then subjected to magnetic bead capture (CpG 
MethylQuest DNA Isolation Kit; EMD MiUipore, BiUerica, Massachusetts, 
USA). Molecules were restricted with Hindlll and ligated into the Hindlll site 
of the pUC19 plasmid. Ligation products were introduced into E. coli strain 
DH5a (ElectroMax; Invitrogen, Carlsbad, California, USA) by electroporation. 
Blue- white selection was used to identify recombinant clones for sequencing on 
an ABI PRISM 377 DNA sequencer (Applied Biosystems, Carlsbad, California, 
USA) using Amersham's DYEnamic ET Terminator Cycle Sequencing Kit 
(Amersham Biosciences, Little Chalfont, Buckinghamshire, United Kingdom). 

One hundred and six clones were sequenced, including 22 from the CA li- 
brary, 22 from the GA library, 20 from the AAC library, and 22 from the ATG 
library. Eighty-two different microsatellite-containing clones were identified 
from the four libraries, and primers were designed for 72 sequences using 
DesignerPCR version 1.03 (Research Genetics Inc., www.lifetechnologies.com) 
and synthesized for 24. These primer pairs were initially tested to verify amplifica- 
tion, determine the optimum annealing temperature, and to establish size ranges 
for later PCR multiplexing, using DNA from eight individuals of P. ellip- 
tica. PCR was carried out in a total volume of 10 \xL, containing approximately 
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10 ng genomic DNA template, Ix PCR Polymerization Buffer containing 
dNTPs (Fisher Biotech, Wembley, Western AustraHa, AustraHa), 0.2 jiM each 
of unlabeled forward and reverse primer (GeneWorks, Hindmarsh, South 
Australia, Australia), 0.5 units Taq DNA polymerase (Fisher Biotech), and 
2 mM of MgCl2 (Fisher Biotech). PCR was carried out in a Veriti 96- Well Thermal 
Cycler (Applied Biosy stems) using the following reaction conditions: an initial 
activation step at 95°C for 15 min; followed by 35 cycles of 95°C for 30 s, 
annealing at 59°C for 90 s, extension at 72°C for 90 s; followed by a final exten- 
sion at 72°C for 15 min. PCR products were separated on a 2% agarose gel 
stained with SYBR Safe DNA Gel Stain (Invitrogen), and fragment sizes deter- 
mined by comparison to a Low DNA Mass Ladder (Invitrogen). 

The nine loci that amplified successfully were then screened on eight indi- 
viduals of P. elliptica to test for polymorphism and trial multiplexing groups. 
Multiplexing was performed on four groups of primers (Table 1) using the 
QIAGEN Multiplex Kit (QIAGEN, Valencia, California, USA) in 12.5-|iL reac- 
tion volumes containing 5-30 ng DNA template, 6.25 jiL QIAGEN Multiplex 
PCR Master Mix, 1.25 |iL Q solution, 0.1 |iM of each forward primer (labeled; 
unique to primer), 0.1 jiM of each reverse primer (unlabeled), and sterile H2O 
to 12.5 jiL. The multiplex PCR was conducted in a Veriti 96-Well Thermocycler 
(Applied Biosy stems) with the following conditions: an initial activation step at 
95°C for 15 min; followed by 35 cycles of 95°C for 30 s, annealing at 59°C for 
90 s, extension at 72°C for 90 s; followed by a final extension at 72°C for 
15 min. Following PCR, samples were diluted 1 : 30 in sample loading solution 
(Beckman Coulter, Brea, California, USA) with the addition of 0.4 jiL fluores- 
cently labeled 400-bp size standard per sample (Beckman Coulter) for capillary 
electrophoresis on a CEQ 8800 Genetic Analysis System (Beckman Coulter). 
Fragment peaks were visualized using CEQ Genetic Analysis System software 
(Beckman Coulter), and fragment (allele) sizes were scored manually. Eight 
individuals of P. longifolia were also genotyped using the above protocol, to test 
for cross-species amplification for these nine loci. 

Genomic DNA of P. longifolia was sent to the Australian Genomic Research 
Facility (AGRF), Adelaide, Australia, for shotgun sequencing on a Titanium 
GS-FLX (454 Life Sciences, a Roche Company, Branford, Connecticut, USA), 
following Gardner et al. (2011). The sample occupied 6.25% of a plate and 
produced 108,806 individual sequences, with an average read length of 367 bp. 
The average GC content of these data was 37.82%. The program QDD version 
1 (Meglecz et al., 2010) was used to screen the raw sequences for eight or 
greater di-, tri-, tetra-, or pentabase repeats, remove redundant sequences, and 
design primers using Primer3 (Rozen and Skaletsky, 2000). Software running 
parameters were set to default values, except PCR product lengths, which were 
set to 80-480 bp. We identified 14.2% of all reads containing microsatellite 



loci. Dinucleotide motifs were the most frequent (79,604), followed by tri-, 
tetra-, and then pentanucleotide motifs (14,550, 1072, and 445, respectively) 
(Meglecz et al., 2012). Primer pairs were designed for 108 loci and, from these, 
we excluded all loci that contained imperfect repeats, had a >2°C difference 
between the forward and reverse primer annealing temperature, short repeat 
motifs within the flanking region or primer sequence, or had poly A/T runs of 
>7 bp as there is an association between a high degree of poly A/T and instabil- 
ity (Li et al., 2002). 

From the remainder, we arbitrarily selected 28 loci and then followed guide- 
lines from Gardner et al. (201 1) for further development; the loci were initially 
trialed for amplification. Initial PCR amplification was carried using the same 
method as for cloned microsatellites, mentioned above. 

The 12 loci that amplified successfully were then screened on eight indi- 
viduals of P. longifolia to test for polymorphism and trial multiplexing groups. 
Multiplexing was performed on four groups of primers (Table 1) using the 
QIAGEN Multiplex Kit in 12.5-jiL reaction volumes containing 5-30 ng DNA 
template, 6.25 |iL QIAGEN Multiplex PCR Master Mix, 1.25 |iL Q solution, 
0.1 |iM of each forward primer (labeled; unique to primer), 0.1 jiM of each 
reverse primer (unlabeled), and sterile H2O to 12.5 jiL. The multiplex PCR was 
conducted in a Veriti 96-Well Thermocycler (Applied Biosystems) with an ini- 
tial activation step at 95°C for 5 min; followed by nine cycles of 95°C for 30 s, 
a 1°C touchdown starting at 65°C for 180 s, and 72°C for 15 s; followed by 25 
cycles at 95°C for 30 s, 56°C for 180 s, and 72°C for 15 s; and then a final exten- 
sion at 60°C for 30 min. Capillary electrophoresis and fragment scoring was 
carried out using the same method as for cloned microsateUites, mentioned 
above. Eight individuals of P. elliptica were also genotyped using the above 
protocol, to test for cross-species amplification for these 12 loci. 

One of the nine loci isolated from microsatellite-enriched libraries for P. el- 
liptica did not amplify rehably in the target species (Table 1). Six out of 12 loci 
isolated from 454 pyrosequencing for P. longifolia were monomorphic, and three 
did not amplify rehably (Table 1). Cross-species ampHfication from P. elliptica 
to P. longifolia was successful for three loci (Table 1). Cross-species amplifica- 
tion from P. longifolia to P. elliptica was successful at one locus (Table 1). 

More than 15 individuals from the two populations of both P. elliptica and 
P. longifolia were then genotyped using nine and six loci, respectively; all loci 
were polymorphic and amplified reliably from the previous screening. Genetic 
diversity parameters and deviation from Hardy-Weinberg equilibrium (HWE) 
were calculated using GenAlEx version 6.4 (Peakall and Smouse, 2006) (Table 2). 
We used MICRO-CHECKER 2.2.3 (van Oosterhout et al., 2004) to check each 
locus for evidence of null alleles, scoring error due to stuttering, and large allele 
dropout. 



Table 1. Characteristics of microsatellite primers developed in Persoonia elliptica and P. longifolia. 















WellRED dye label 


GenBank 


Locus 


Method^ 




Primer sequences (5 '-30 


Repeat motif 


Allele size range (bp)^ 


(Multiplex marker set) 


accession no. 


PeA102 


GIS 


F 
R 


ACCCACCAGCTAACTCTTTTG 
CCATCTGATCGTACCAACGTA 


(CA)i4 


94-138 


D3 (3) 


JX987704 


PeA103 


GIS 


F 
R 


CCTAATAACCCCACCCTACTAC 
AGGTCCTTTCACACTTTTGAT 


(GT)i5 


243-299; 243-289 


D3 (2) 


JX987705 


PeA106 


GIS 


F 
R 


TCCAACTGGATGCCAACTC 
TCTTGGGATCAGCGTATGC 


(GC)5(AC)i5 


179-187; 167-169 


D2(l) 


JX987706 


PeB103 


GIS 


F 
R 


AATGGGGTTGCTGTGTCT 
TCCCTTCTCTTCCACCAC 


(GA)i9 


272-314; 270-282 


D4 (2) 


JX987707 


PeB106 


GIS 


F 
R 


TGCTTCTTGTATGCCAAAGAG 
CCTTCTAAGGGAGATGCAGAG 


(GTT)9 


164-226 


D2 (3) 


JX987708 


PeB2 


GIS 


F 
R 


TTCGCTTCCTTCTTTGTTAAAG 
AGGGAGGTATCACATGAGGAC 


(CT)i9 


285-313 


D4(l) 


JX987709 


PeC104 


GIS 


F 
R 


C T T GAGC AT GAGC T T GAGAT C 
AT C AAC C AT GAGGAGAAGAGAG 


(GTT)9 


146-168 


D4(l) 


JX987710 


PeC2 


GIS 


F 
R 


GATCCTTAGGTTTGGTTTCTTG 
TGCCCTATCAAGCATAGTTATG 


(GTT)io 


235-267 


D3(l) 


JX987711 


P113 


454 


F 
R 


GAGGCCATATAGGGTGGTGA 
CACAGGATCATAGTTGCCTCG 


(AT)8 


na; 110-132 


D2 (5) 


JX987701 


P115 


454 


F 
R 


ATAATGGGCCTAATGATGCG 
ACTCCCTAGGCTATGCCTCC 


(GA)9 


145-163; 135-181 


D3 (3) 


JX987702 


P118 


454 


F 
R 


TGGACACTTCAGGAGCTAACA 
GGGAGGTGTATCCACTCCTTG 


(GA)i5 


na; 188-206 


D2 (5) 


JX987703 



^The two development methods were (a) microsatellite-enriched library construction by Genetic Identification Services (GIS; La Canada Flintridge, 
California, USA; http://www.genetic-id-services.com/) for P. elliptica and (b) 454 GS-FLX shotgun sequencing for R longifolia. 

^"na" indicates that the primer was not used for P. elliptica; for primers that are used in both species, values for P. elliptica are listed first. 
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Table 2. Results of primer screening in two populations each of Persoonia elliptica and P. longifolia.^ 



Locus 








P. ellipt. 


ica 














P. longifolia 










Pe-AW = 15) 






Pe-AE = 15) 






Pl-MD ( 


n - 30) 






Pl-CE ( 


n = 25) 




A 


"o 


ZJ 




A 


ZJ 

^0 


ZJ 

"e 


HWE^ 


A 


Ho 




rl W Jj/ 


A 

/± 


Ho 


H, 


HWEb 


PeA102 


10 


0.46 


0.84 




10 


0.73 


0.82 


ns 


















PeA103 


7 


0.80 


0.78 




5 


0.80 


0.65 


ns 


6 


0.70 


0.68 


ns 


6 


0.76 


0.68 


ns 


PeA106 


3 


0.47 


0.42 


ns 


4 


0.67 


0.61 


ns 


3 


0.07 


0.07 


ns 


2 


0.04 


0.04 


ns 


PeB103 


14 


0.87 


0.86 


ns 


10 


0.67 


0.80 


ns 


6 


0.40 


0.46 




5 


0.60 


0.68 


ns 


PeB106 


13 


0.69 


0.84 




11 


0.87 


0.86 


ns 


















PeB2 


8 


0.93 


0.84 


ns 


10 


0.93 


0.88 


ns 


















PeC104 


4 


0.73 


0.66 


ns 


4 


0.87 


0.68 


ns 


















PeC2 


6 


0.93 


0.78 




8 


0.80 


0.81 


ns 


















P113 


















3 


0.30 


0.32 


ns 


5 


0.56 


0.68 




P115 


10 


0.86 


0.85 


ns 


8 


0.93 


0.82 


ns 


13 


0.73 


0.84 


ns 


10 


0.76 


0.79 




P118 


















6 


0.57 


0.68 


ns 


6 


0.88 


0.68 


ns 



Note: — = primer was not genotyped; A = number of alleles; = expected heterozygosity; = observed heterozygosity; HWE = Hardy-Weinberg 
equilibrium; n = number of individuals sampled. 

^All values are based on samples from populations in southwestern Austraha (Pe-AE: 31.572°S 116.181°E, Pe-AW: 31.613°S 116.151°E, Pl-MD: 
32.176°S 116.245°E, Pl-CE: 33.356°S 116.345°E). 

''Significant departure from HWE is shown as: *P < 0.05, **P < 0.01, *** P < 0.001; ns = not significant. 



In P. elliptica, between three and 14 alleles per locus were found. Observed 
(//J and expected (H^) heterozygosities ranged from 0.46 to 0.93 and 0.42 to 
0.88, respectively (Table 2). For P. longifolia, between two and 13 alleles per 
locus were found. and ranged from 0.04 to 0.88 and 0.04 to 0.84, respec- 
tively (Table 2). Significant departure from HWE was detected in four of the 
nine loci for P. elliptica and three of the six loci for P. longifolia (Table 2). No 
loci showed significant null allele frequencies, large allele dropout, or evidence 
of scoring error due to stuttering. 

CONCLUSIONS 

The microsatellite markers developed here will enable an ex- 
amination of population genetic patterns and realized dispersal of 
pollen and seed in R elliptica and P. longifolia. The effort re- 
quired in developing microsatellite markers has long been an 
important consideration, especially for low-budget studies on 
nonmodel organisms (Squirrell et al., 2003). Next-generation se- 
quencing techniques such as 454 pyrosequencing have now been 
widely used for the development of microsatellite markers more 
rapidly and cheaply than traditional cloning approaches (Gardner 
et al., 2011; Malausa et al., 2011). Our experience suggests that, 
although there were efficiencies to be gained from employing 
next-generation sequencing for microsatellite marker develop- 
ment, particularly in the identification of microsatellite loci and 
primer sequences, it is not a panacea. A significant investment of 
time and effort in the screening, as well as optimization of poly- 
morphic microsatellite markers, is still required with this rela- 
tively new approach to marker development. 
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